Geometric Landau-Zener interferometry 
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We propose new type of interferometry, based on geometric phases accumulated by a periodically 
driven two-level system undergoing multiple Landau-Zener transitions. As a specific example, we 
study its implementation in a superconducting charge pump. We find that interference patterns 
appear as a function of the pumping frequency and the phase bias, and clearly manifest themselves 
in the pumped charge. We also show that the effects described should persist in the presence of 
realistic decoherence. 



A driven quantum two-level system traversing an 
avoided energy-level crossing can undergo nonadiabatic 
transitions, known as the Landau-Zener effect. If more 
than one crossing is involved and the dynamics is overall 
coherent, then transition paths can interfere according 
to the different phase accumulated by the ground and 
excited-state wavefunctions between subsequent cross- 
ings. This phenomenon, sometimes referred to as 
Landau-Zener-Stiickelberg (LZS) interferometry pQ, was 
first observed in atomic and optical systems, and recently 
proposed [2] and measured also in superconducting qubit 
systems [3HZ]- In all these realizations, the system is 
driven in such a way that the interference effects have 
a purely dynamical nature. In general, though, a quan- 
tum state subject to steered evolution acquires both a 
dynamic and a geometric phase. While the study of ge- 
ometric phases in solid-state systems is an active field of 
research [5MlO|. their relevance to LZS interferometry has 
so far been unexplored [TTj . 

In this Letter, we elucidate the link between LZS inter- 
ference and geometric phases, opening new possibilities 
for the geometric control of quantum systems. Our re- 
sults apply to a broad range of devices, namely those for 
which the parametric driving possesses a nontrivial geo- 
metric structure and the induced energy-gap modulation 
presents multiple avoided crossings. As a pertinent ex- 
ample, we consider a superconducting charge pump, the 
Cooper-pair sluice [T^] . The connection between Cooper- 
pair pumping and geometric phases was highlighted in 
previous theoretical works both for the abelian [T3J H3] 
and nonabelian |15j case, yet always in the adiabatic 
limit, where the system stays in the instantaneous ground 
state and excitations are treated as small corrections [TB] ■ 
We instead consider higher frequency regimes and predict 
the appearence of interference patterns depending on the 
pumping frequency and the superconducting phase bias, 
the latter embodying the geometric contribution to in- 
terference. We then show that LZS resonances directly 
manifest themselves in the pumped charge, which is an 
advantage of using a charge pump rather than a conven- 
tional qubit as an interferometer. Finally, we introduce 
decoherence in our model and show that interference ef- 
fects are still detectable. This should make our proposal 
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FIG. 1. (Color online) (a) Schematic drawing of the "sluice", 
(b) Effective magnetic field corresponding to the pumping cy- 
cle considered in this Letter, (c) Adiabatic (instantaneous) 
energy versus time. Avoided level crossings occur at times 
ti = T/4 and ti = 3T/4. Green and red arrows outline two 
possibly interfering paths. 



feasible for experimental observation. 

The Cooper-pair sluice, schematically shown in 
Fig. 1(a), consists of a superconducting island coupled 
to the leads by two superconducting quantum interfer- 
ence devices (SQUIDs), whose Josephson energies Jl,r, 
can be tuned by changing the magnetic fluxes $>l,r- A 
gate electrode capacitively coupled to the island is used 
to induce a polarization charge n g on the latter, thereby 
providing a third control parameter. During a pumping 
cycle, the parameters are steered so as to couple the is- 
land to the left lead, attract a Cooper-pair, switch the 
coupling to the right lead, and release the Cooper-pair 
[T7] . We will assume that the superconducting phase dif- 
ference <j> across the device is kept constant. This can 
be achieved by shunting the sluice with a large Joseph- 
son junction. In this case, the switching statistics of the 
additional junction also provides a way to measure the 
pumped charge [18] , 

The device is operated in the Coulomb-blockade 
regime E c > J max , where J max = max{ J L ,J R } and 
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Ec is the charging energy of the island. The system dy- 
namics is then best described in the basis of eigenstates 
of charge on the island. Also, as long as n g stays close 
to the degeneracy point 1/2, only two such states are 
relevant, namely those with zero and one excess Cooper 
pair on the island. This allows us to use a pseudo-spin 
formalism and write the sluice Hamiltonian as H = a- B, 
where {<Ji\ are the Pauli matrices and the effective mag- 
netic field B has components 
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B x {t) = J+(t) cos 

B y (t) = J- (i) sin I , 
B z {t) = E c [1/2 - n g (t)] , 



(1) 

(2) 
(3) 



where we put J±(t) = Ji(i)± Jn(t). As B is steered along 
the path shown in Fig. 1(b), it spans a solid angle which 
is responsible for the geometric effects under discussion. 
This situation is clearly different from that considered in 
e.g. Refs. [3JH], where B moves on a definite plane (say, 
x-z), leaving no room for non-trivial geometric effects to 
take place. 

In Fig. 1(c) we plot the energies of the adiabatic states 
\g), |e) as a function of time for a pumping cycle, ob- 
tained by instantaneous diagonalization of the Hamilto- 
nian. The avoided level crossings at t — T/4, 3T/4 (T 
is the pumping period) correspond to the gate charge 
crossing the degeneracy point. The probability of a nona- 
diabatic (Landau-Zencr) transition at such a crossing is 
given by Plz = e~ 2 ^ s , where the adiabatic parameter 8 
depends on the velocity at which the crossing is traversed 
and on the energy gap at the crossing [T . For our case, 
8 = irJmaxI (48Ecfighv), where n g = max{n ff } — 1/2 and 
v = l/T. 

In the limits Ec S> J max and hv < Ecn g , nonadia- 
batic transitions are strongly localized at level crossings. 
The system dynamics can thus be seen as a sequence of 
adiabatic evolutions and localized transitions. For this 
reason, the calculation is most conveniently performed in 
the adiabatic basis {\g), |e)}. In the so-called adiabatic- 
impulsc model (TJIHJ], Landau-Zencr tunneling at anti- 
crossings is treated as instantaneous and described in the 
adiabatic basis by a transfer matrix of the form: 



N LZ = 
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where (ps = 5 (log 5— l) + argr(l — iS) — ir/4 is the impul- 
sive phase acquired by the adiabatic states in traversing 
the crossing (T is the gamma function) |U |5U] . 

For each adiabatic segment j — 1,2,3 in Fig. |ljc) , 
evolution from time to tj is described by a diagonal 
matrix of the form: Uj — exp [i(pj<r z ] , where ipj is the 
total phase difference acquired by the adiabatic states. 
The latter can be written as ipj = £j + jj, where we 
have distinguished a dynamic and a geometric (7^) 



contribution. The dynamic phase difference £j is given 
by: 
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dtV\H 11 -H 22 \ 2 +4\H u 



(5) 



By contrast, the gauge-invariant, noncyclic geometric 
phase difference jj can be calculated as [2Tj : 



dt 



(6) 



This phase is uniquely determined by the path drawn by 
the system in parameter space, and reduces to the Berry 
phase for cyclic ground-state evolution [2"2"] . 

Putting things together, the evolution operator over a 
period can be calculated as 



U = U Z N LZ U 2 N LZ XJ 1 = U 3 



-0* 



Ui , (7) 



where 



a = [(1 - Plz) 



-P LZ e 



-i£2-*72l 



P = 2y/P LZ (l-P LZ ) COS (cp S + 6 + 72) 



(8) 
(9) 



U\ and U3 play no role in the upcoming resonance con- 
dition and will not be considered further. From ^7b, we 
can calculate the excitation probability after one period 
starting from the ground state. This is given by: 

P = 4P lz (1-P lz )cos 2 (<as + 6 + 72) ■ (10) 

This probability oscillates between and APl Z (\ — Plz) 
as a function of the accumulated phase. In the fast- 
passage limit, 8 <^ 1, we can approximate ips ~ — 7r/4 
(in the adiabatic limit 8 — > 00 and (ps — > —ir/2). 

We now make our discussion specific by considering 
the pumping cycle of Fig. 1(b) |23]. Up to the first order 
in [J ma x/(E c n g )} 2 , we find: 
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As expected, the dynamic phase £2 is inversely pro- 
portional to the pumping frequency v. By contrast, the 
geometric contribution 72 does not depend on i>, and 
in this particular case equals half the superconducting 
phase bias </>. We have thus derived a resonance condi- 
tion involving the superconducting phase bias <j) and the 
pumping frequency v. In particular, in the region where 
(ps attains a constant value, the resonances drift in the 
<p — v plane as branches of hyperbolae. 

This analysis predicts the position of resonances and 
explains their origin. Its regime of validity lies in be- 
tween the strictly adiabatic and the fully nonadiabatic 
one. As a matter of fact, a lower bound for the pumping 
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FIG. 2. (Color online) (a) Pumped charge (in units of e) after 
one cycle versus phase bias and frequency. The parameters 
are: E c /k B = 2.5 K, J max = O.lEc, 0.3 < n g < 0.7. Dashed 
lines enclose the regions where the ground state population 
at the end of the cycle is at least 0.9. Dotted lines have the 
same meaning but they are calculated according to ( 10 1 and 
in the fast-passage approximation ips ~ — tt/4. (b) Ground- 
state population versus time for a case of destructive (dashed 
line) and constructive (solid line) interference. The pumping 
frequency is 1.56 GHz for both cases, the phase biases are 
0.22 and 3.36, respectively. 



frequency is set by the requirement for time evolution to 
be coherent over one pumping cycle. On the other hand, 
at frequencies comparable to the adiabatic level spacing 
(hv w Ecfig) transitions are no longer restricted to the 
degeneracy points and the adiabatic- impulse model is ex- 
pected to break down. 

The superconducting phase bias <p enters the resonance 
condition through the geometric phase accumulated be- 
tween subsequent transitions. This relationship is trivial 
for the case considered, as the geometric phase is simply 
proportional to <j>. Yet, this example clearly illustrates 
the role of geometric phases in Landau-Zcncr interfer- 
ence. In particular, by choosing the pumping frequency 
so that (ps + £2 is an integer multiple of w, the dynamic 



contribution in (10 1 is washed out, resulting in a purely 
geometric Landau-Zener interference effect. 

We now proceed to show that the predicted resonances 
manifest themselves in the charge pumped by the device, 
thus providing the most straightforward way of observing 
them. To do so, we first obtain the full system dynam- 
ics from numerical solution of the Schrodinger equation. 



We then calculate the pumped charge by integrating the 
instantaneous current operator [24 . In Fig. 2 (a), we 
plot the pumped charge over a period versus the phase 
bias <p and the pumping frequency v. The parameters 
are chosen so as to be consistent with our model. In 
particular, microscopic excitations in the superconduct- 
ing circuit can be neglected provided hf e ff -C A, where 
f e ff is the effective frequency of the driving fields and 
A the superconducting gap. Furthermore, for small val- 
ues of fig the system is sufficiently anharmonic for the 
two-level approximation to hold in the given frequency 
range. The lines drawn on top of the image plot cor- 
respond to 90% probability of the system being in the 
ground state at the end of the cycle. Dashed lines are cal- 
culated numerically, dotted lines according to (fl0|. The 



strong correlation between the ground state population 
and the pumped charge demonstrates the possibility to 
access interference patterns simply by measuring the lat- 
ter. Moreover, the accuracy of the approximations made 



in deriving (10 1 is confirmed by the good agreement be- 
tween analytical and numerical calculations. 

In Fig. 2 (b), we show the time evolution of ground- 
state populations for one case of constructive and one of 
destructive interference. In both cases there is a popula- 
tion transfer to the excited state after the first crossing. 
Yet, while constructive interference (solid line) enhances 
the excitation after the second crossing, destructive in- 
terference (dashed line) brings the system back to the 
ground state. In particular, this implies that for a given 
pumping frequency, the phase bias can be chosen so as 
to pump a significant fraction of Cooper pairs (about 0.5 
in this case) even in the nonadiabatic regime. 

A complementary and instructive way to understand 
these features is provided by Floquet analysis [25] • In 
fact, we can explicitly calculate the quasienergy spectrum 
by diagonalizing the evolution operator U in (j7J| . We find 
that destructive resonances occur at exact quasienergy 
crossings, where time evolution over a period is trivial 
and tunneling between adiabatic states is dynamically 
frozen. This phenomenon is known as coherent destruc- 
tion of tunneling [26j . At the opposite end, constructive 
interference enhances such transitions, resulting in Flo- 
quet states being the maximal mix of the adiabatic ones. 
This is revealed in the quasienergies as the opening of a 
gap, similarly to a time-independent system with a cou- 
pling interaction switched on. 

The LZS interferometry discussed above is a unitary 
and coherent process. However, in any experimental im- 
plementation the undesired and unavoidable coupling to 
external degrees of freedom leads to decoherence and thus 
it affects the observed pumped charge. This fact pre- 
cludes a measurement of coherent LZS pumping over a 
great number of cycles and must be taken into account 
in any experimental proposal. To this end, we will now 
discuss the case in which the system is affected by charge 
noise. The dissipative dynamics of the system is numer- 
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ically obtained from the master equation including the 
driving field and the environment. The latter is described 
by a bath of harmonic oscillators with ohmic spectrum 
at zero temperature [27J US] ■ To be able to detect the 
effect of coherence loss, the pumping period must be 
smaller than the expected decoherence time. By choos- 
ing a superconducting island with high charging energy 
Ec/ks = 5 K, this condition is fulfilled at frequencies as 
low as 0.5 GHz, which is still in the regime of validity of 
the master equation [271 129] . 

In Fig. [3] (a) we present the expected pumped charge 
in the absence of noise over 50 consecutive pumping cy- 
cles for constructive (squares) and destructive (circles) 
interference. The results are readily interpreted after 
Fig. [2] (b) : When interference is destructive, the system 
starts every cycle in the ground state, so that the pumped 
charge is constant. On the contrary, constructive inter- 
ference allows the system to be in different superpositions 
of the ground and the excited state, and this is reflected 
in an oscillating pumped charge from cycle to cycle. 

When we include the environmental effects as in 
Fig.[3jh), the behavior of the pumped charge in the con- 
structive case changes dramatically, as oscillations are 
quickly damped due to the loss of coherence. At the 
same time, the pumped charge in the destructive case 
is only slightly affected by the environment. This is a 
direct result of the fact that the system stays mainly in 
the ground state, which was shown to be robust against 
decoherence [27J |2S]. This suggests that the detrimen- 
tal effects of decoherence can partly be overcome by a 
suitable choice of the phase bias. Most importantly, we 
notice that the large difference in the pumped charge 
makes it possible to clearly distinguish between the two 
cases. 

In conclusion, we have proposed new type of Landau- 
Zener interferometry, based on geometric phases. Specif- 
ically, we have demonstrated this technique in a super- 
conducting charge pump. For that case, we have shown 
that the geometric effects are controlled by the super- 
conducting phase bias across the pump. They can be 
detected by measuring the pumped charge, and should 
persist in the presence of realistic decoherence. 
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